STRUCTURAL INVARIANCE AND THE ENERGY SPECTRUM 
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We extend the application of the concept of structural invariance to bounded time independent 
systems. This concept, previously introduced by two of us to argue that the connection between 
random matrix theory and quantum systems with a chaotic classical counterpart is in fact largely 
exact in the semiclassical limit is extended to the energy spectra of bounded time independent 
systems. We proceed by showing that the results obtained previously for the quasi-energies and 
eigenphases of the S-matrix can be extended to the eigenphases of the quantum Poincare map 
which is unitary in the semiclassical limit. We then show that its eigenphases in the chaotic case 
move rather stiffly around the unit circle and thus their local statistical fluctuations transfer to 
the energy spectrum via Bogomolny's prescription. We verify our results by studying numerically 
the properties of the eigenphases of the quantum Poincare map for billiards by using the boundary 
integral method. 



I. INTRODUCTION 

In previous papers two of us developed the concept of 
structural invariance for periodically driven |l| as well as 
for scattering systems [HJ3| in order to establish the con- 
nection between random matrix theory on the one hand 
and quantum systems with a chaotic classical counter- 
part on the other. Basically we recovered for such sys- 
tems the result predicted for the long-range behaviour 
of the spectral two-point function in a seminal paper by 
Berry [H and demonstrated numerically and experimen- 
tally in a large, and growing, number of examples |0,E3[: 
Namely, the fluctuation properties of quantum systems 
whose classical counterpart is chaotic behave as those 
of an appropiately chosen ensemble of random matrices. 
There have been numerous attempts to account for this 
connection. Among the most notable are periodic orbit 
theory, which shows how certain random matrix proper- 
ties arise from a semiclassical formula connecting the den- 
sity of states with a given sum over periodic orbits. An- 
other approach, which has achieved some remarkable suc- 
cesses recently, is the supersymmetric approach, which 
has been able to show how such phenomena arise in the 
context of disordered systems g. For chaotic systems 
without disorder, however, the results obtained in this 
manner are still far from convincing nM- The approach 
taken by structural invariance is in some sense different, 
since it takes an explicitly probabilistic point of view, as- 
serting only that a chaotic system has RMT properties 
with probability one and only speaking in terms of en- 
sembles of systems. [|]|lj,^. This approach has already 
led to unexpected predictions, which since have been con- 
firmed numerically ||9|] and later by arguments taken from 
periodic orbit theory [ p^ . 

A weak point of our approach was the fact that its 
application was in principle restricted to the discussion 
of eigenphase statistics. In this paper, we shall use the 
semiclassical quantization of the Poincare map and its 
relation to the energy spectrum in order to transfer the 
statistics of eigenphases to that of eigenvalues. The fol- 
lowing relation, due to Bogomolny EG] 



det(l - T{En)) = 



(1) 



shows that we obtain an energy eigenvalue En each time 
one of the eigenphases of the matrix T (operator asso- 
ciated to the quantum Poincare map) takes the value 
unity; thus a stiff rotation of the unit circle on which 
the eigenphases are located would transfer their spectral 
fluctuations to those of the energy eigenvalues. In this 
paper we shall show that the assumption of fairly uni- 
form speed of the eigenphases as a function of energy is 
an excellent approximation for chaotic Poincare maps. 

The structure of this paper is as follows: First, in Sec- 
tion 2, we give a short account of structural invariance 
in the most general terms possible. In Section 3, we 
review how the connection between classical behaviour 
and quantum statistics is effected in the case of canon- 
ical maps. In Section 4, we proceed to generalize the 
analysis to Poincare maps. There we show, using Bo- 
gomolny's semiclassical quantization, that the spectral 
statistics of the eigenphases of the Poincare map essen- 
tially carry over to those of the eigenvalues. In Section 
5, we illustrate the results by means of a numerical com- 
putation of the stadium billiard. 



II. STRUCTURAL INVARIANCE 

The point we want to discuss in this Section concerns 
the construction of a reasonable ensemble given an in- 
dividual system. Such problems are not really new to 
physics: indeed, whenever one attempts to describe an 
individual system by an ensemble, as occurs for example 
in statistical mechanics, the problem of deriving the en- 
semble from the individual system arises. In such cases 
it is well known that one first needs to identify a set of 
relevant properties. Once this has been done, we define 
the structural invariance group of the object as the group 
of all transformations which transform the given object 
into one with the same relevant properties. 

First, let us present a general method to derive such 
an ensemble. Let us consider a set of objects with given 
"relevant" properties — the choice of these will turn out to 
be crucial, though it is usually dictated from the prob- 
lem at hand — and a specific object A taken from this 
set. Furthermore, let the various objects be transformed 



among each other through the action of a group Q. The 
object A has a certain number of the "relevant" proper- 
ties. We define 7i, the structural invariance group of A, 
in the following manner: Let Ti. be that subgroup of Q 
which leaves all properties of A invariant. If we now act 
on A with all the transformations of 7i, we generate a 
collection n of objects, which has as additional structure 
that of a homogeneous space under the group 7i. li Ti. 
additionally has an invariant measure, then a measure on 
f2 can be generated in a natural way as follows: To each 
subset E of 17 one can associate the set of transformations 
T in Ti. which map ^ in S. The Haar measure of T can 
then be taken as the measure of S. This is independent 
of the choice of A, because of the invariance of the Haar 
measure. Furthermore, it is a measure on fl that is in- 
variant under the group action and it is therefore clearly 
singled out. 

Let us first consider a trivial example: Consider the set 
of all infinite binary sequences of ones and minus ones. As 
relevant properties, we take the value of all finite range 
correlation functions. By this we mean the following: 
Consider the sequence (s^). The two-point correlation is 
then given by 



Ci = lim — — 
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and more general correlations defined through the ob- 
vious generalization. As transformations acting on bi- 
nary sequences, we consider all possible transformations. 
However, the resulting group is highly pathological and 
has presumably no Haar measure. One way to proceed is 
to limit ourselves to large but finite sequences. The group 
Q then consists of all transformations mapping arbitrary 
finite sequences on other arbitrary finite sequences. If the 
correlations of the original sequence are sufficiently small, 
then this property is respected by almost all transforma- 
tions of G, that is all but a set of very small measure. 
The structural invariance group Ti is then the full group. 
From this follows that we can consider the original se- 
quence as a representative element of the full ensemble 
of all binary sequences, in which all sequences are taken 
to be equally probable. On the other hand, if there is a 
clear predominance of ones, say, but no further correla- 
tions, then we can take Ti to be the set of all maps which 
permute the elements of the original sequence. The set il 
is then the set of all sequences which arise from the origi- 
nal one through permutation, all sequences being equally 
probable. Again, it is reasonable to assume that the orig- 
inal sequence is indeed a representative member of the 
ensemble. In this sense, whatever property holds with 
probability one in the ensemble that has been constructed 
in this fashion will also presumably hold for the original 
element. Note that this approach does not always work: 
For example, if I find a non-zero two-point correlation, 
it is at least not evident how to define a group Ti that 
respects this property. A better approach might then be 
to define a non-invariant measure on Q that takes this 
correlation into account and that respects it with proba- 
bility one. However, it will turn out that we shall be able 



to define a structural invariance group in the majority of 
the cases we shall be dealing with. 



III. CANONICAL MAPS: ASOCIATED 
QUANTUM AND CLASSICAL ENSEMBLES 

We now wish to apply the above general approach to 
the particular case of bijective canonical maps on a com- 
pact phase space. To this end, we must define the various 
ingredients involved in the construction. The set of all 
objects is the set of all bijective canonical transformations 
on a compact phase space F. The specific object wc start 
from is a given such transformation Cq. As for the "rel- 
evant" properties, their choice is dictated by the nature 
of the problem. In the foUowoing, we shall principally 
be interested in the short-range statistical properties of 
the eigenphases and eigenvalues of the associated unitary 
and hermitian operators. This corresponds to the long- 
time properties of the classical time evolution. For this 
reason, we shall consider only such properties to be "rele- 
vant" as remain invariant under arbitrary iteration of the 
map. As for the group G, we shall take it to be the direct 
product of the group of all bijective canonical canonical 
transformations with itself, with the following action on 
the set of all bijective canonical transformations: 



PCo,,Cfi '■ C —> C'aCC'fj. 



(3) 



To complete the construction, it would be necessary to 
have a Haar measure on the group Q. As far as we know, 
the existence of such a measure has neither been proved 
nor shown to be impossible. While we believe that such 
a measure does indeed exist in the case of compact phase 
spaces, we shall not make siuch an assumption, but rather 
make the following consideration: After all, what we are 
really interested in are not the (classical) properties of 
the canonical maps, but rather those of the corresponding 
quantum unitary maps. We therefore need to quantize 
both Co and the set of all canonical transformations. To 
this end, we first need to know the quantum equivalent 
of the phase space F. Since F is compact, its quantiza- 
tion will be a finite dimensional Hilbert space, with the 
dimension N going to infinity as ?i — > 0. Note that in 
general not all values of h yield a quantization of F, but 
only such for which cells of size {2TThY fit in an integer 
nuber of times in the phase space F. We now proceed 
to quantize all canonical maps using some quantization 
technique. Two things shoud be noted here: First, there 
is no unique way to quantize a canonical map. Second, 
whichever way is used to associate a unitary map to a 
canonical map, it can never be such that the quantiza- 
tion of the produvct of two canonical maps is the prod- 
uct of the two quantized maps. Such a result can in 
deed be obtained, but only approximately, in the limit 
N -^ oo. After having constructed il in the set of all 
classical maps, we can proceed to define fig to consist 
of the set of all unitary maps corresponding to a given 
quantization of the maps in il. In a similar way, we can 
translate the group Q in the quantum domain, where it 



becomes simply U{N) x U{N), where N is the dimen- 
sion of the Hilbert space. The finiteness of TV is a crucial 
point here, and this is precisely the point at which we 
use the fact that the phase space F is compact. From 
this follows that there exists a (unique) Haar measure on 
the quantized version of Q and hence, for every closed 
subgroup 7i as well. 

There is an apparent problem with the above program, 
however: Since the "translation" from classical to quan- 
tum language is by no means unique, one cannot assign to 
each classical canonical map a well-defined unitary map. 
Nevertheless, since we are only interested in the semi- 
classical limit, this problem is not really severe, since all 
different possible choices will be very close to each other. 
One might worry, further, that an approximate transla- 
tion of the group Ti. from the classical to the quantum 
domain might lead one to lose the group property. But 
this is not the case, since the group Ti. is defined through 
the invariance of the properties of Co. This definition 
can itself be carried over to the quantum domain and the 
group property is then trivially maintained. 

Another point that can be raised is the following: 
Most canonical maps of interest are "simple" or smooth, 
whereas there is no reason to believe that such smooth- 
ness is characteristic of "arbitrary" canonical maps, if 
this can be given a meaning. Certainly, for unitary maps, 
it is easy to see that arbitrary unitary maps in a (high- 
dimensional) Hilbert space representing the semiclassical 
counterpart of a given classical system with probability 
one fail to correspond to any reasonable classical canoni- 
cal transformation. In the light of these facts, one might 
worry that our reasoning only applies to such systems as 
are already arbitrarily complex at the classical level, and 
are therefore inapplicable in the cases of interest, where 
the map is smooth. We argue that this is incorrect, for 
the folowing reason: Since we are mainly interested in 
short-range spectral statistics and since these are domi- 
nated by the high iterates of Co , we are therefore in fact 
more interested in C^ for iV ^ 1 than in Cq itself. But 
C^ is indeed highly irregular, since Co is chaotic. For 
this reason, it is presumably qualitatively not different 
from C^, where C is an "arbitrary" canonical transfor- 
mation. While such a claim is obviously hard to substan- 
tiate, particularly in the absence of a measure on the set 
of canonical transformations, which would give a precise 
meaning to the word "arbitrary" in the above phrase, it 
is very reasonable and intuitively appealing. 

A related issue arises with respect to integarble and 
near-integrable systems: Since sufficiently weak smooth 
perturbations of an integrable system has invariant tori 
which cover a non- vanishing measure of phase space, it 
might appear that the set of all such systems should have 
positive measure in any ordinary sense of the word, in 
contradiction to our claim, which is that almost all sys- 
tems are fully chaotic. This issue is resolved by noting 
that a generic perturbation in our sense of the word will, 
in general, not staisfy the smoothness requirements of 
the KAM theorem, so that any amount of such a generic 
perturbation will lead to completely irregular dynamics 
for sufficiently high iterates of the map. 



Finally, it remains to see what are the possible special 
properties of Cq. Since we agree that the only relevant 
properties are those that remain invariant under arbi- 
trary iterations of the map, we are left with very few 
possibilities. Here are the three that we are aware of: 
— Time-reversal Invariance: If there exists an anti- 
canonical map T (that is, a map which maps the canoni- 
cal symplectic two-form on its negative) such that T^ = 1 
and such that 



TCoTCo = 1, 



(4) 



then the map Cq is said to be time-reversal invariant 
(TRI). In this case, we take H to consist of the pairs 
{C,TC^^T). These form a group, and it is straightfor- 
ward to calculate that it leaves the TRI property invari- 
ant. As such, it is the structural invariance group of an 
arbitrary Cq which has no specific property apart from 
TRI. From this follows that il is simply the set of all TRI 
bijective canonical maps, as was in fact to be expected. 
On the other hand, we are not (yet!) able to define a 
measure on Q., since neither Q nor Ti have known Haar 
measures. 

This can be readily translated into quantum-mechanical 
language if we note that in the absence of spin, T can be 
represented as complex conjugation in quantum mechan- 
ics, so that H is given by the subgroup consisting of the 
pairs of the type {U, C/*), where C/* is the transpose of U. 
— Discrete Symmetries: Let Co commute with a group 
C of canonical transformations Q. This means that the 
evolution given by Co has symmetries, given by the maps 
Q. This property is indeed preserved for arbitrary iter- 
ations of Cq. The structural invariance group Ti. is then 
given by the subgroup of all pairs (C, C) such that both 
C and C commute with all elements of G. Under these 
circumstances, the action of Ti on Cq respects the prop- 
erty of commuting with C and generates as a set il the 
set of all bijective canonical maps that commute with C. 

Translating into quantum mechanics, this means that 
the quantum version of Ti consists of all pairs of unitary 
matrices that commute with the quantization Uq of the 
maps Q. This means that the unitary matrices acting 
to the left and right of Uc can be put in block diagonal 
form in such a way that each block belongs to a given 
irreducible representation of G. Thus, the ensemble gen- 
erated has a well defined block structure, in which all 
blocks statistically independent. On the other hand, if 
a given block transforms according to a rf-dimensional 
representation of G, the eigenvalues of the corresponding 
block will be d-fold degenerate. 

— Systems with mixed dynamics: A related problem 
arises with systems having invariant tori, though not, as 
we shall see, cantori. Indeed, invariant tori are structures 
which are preserved under arbitrary iterations. Further- 
more, at least in the most frequently studied case of two 
degrees of freedom, they also lead to an absolute division 
of the phase space in disjoint parts. The structural in- 
variance group should therefore respect these character- 
istics. The canonical maps belonging to Ti should there- 
fore have the same invariant tori as the map Co iself. It 



will therefore automatically also leave the same chaotic 
regions invariant. 

Translating this into quantum mechanics, we see that 
to each torus there corresponds a unique wave function 
via the WKB prescription. The above condition there- 
fore means that the group of unitary maps must act on 
all these integrable eigenfunctions through a phase fac- 
tor only. From this follows that the eigenphases in this 
case can be separated in an integrable spectrum, which 
consists of statistically independent eigenphases — since 
they are acted upon by the group U{1) with its Haar 
measure — and a chaotic part, which can be treated in the 
usual way described above, and therefore corresponds to 
a COE or a CUE for each chaotic region. Note that this 
follows from the preceding remark: Indeed, when two 
chaotiic regions are separated by an invariant torus (or 
an integrable region) , the index of the chaotic region is a 
discrete conserved quantity and acts exactly in the same 
way as a discrete symmetry. 

This leads to various questions concering the relevance 
of other phase space structures, such as homoclinic and 
heteroclinic tangles, barriers to transport and cantori. 
That barriers to transport only define separate phase 
space regions up to a given number of iterations is clear. 
For this reason, we should certainly not consider it now, 
since we are limiting ourselves to the simplest, fully semi- 
classical case. However, it is clear that such issues will 
eventually need to be addressed, since structures of this 
type can notoriously cause strong localization effects. 
While it is aggreed that these effects are transient when 
phase space is compact, that is the effects disappear at 
sufficiently small h, they are nevertheless extremely im- 
portant in practice. Similar arguments allow to rule out 
the other structures as well: Since they are all invariant 
sets, it might be argued that they survive under arbi- 
trary iteration. This is true as far as it goes, but they 
are complex sets, with fractal structure, which quantum 
mechanics can never fully capture. If we take a finite 
approximation to such an invariant set, however, it will 
become ever more complex under iteration and there- 
fore, from the point of view of quantum mechanics, it 
will eventually become irrelevant. 

— Non-primitive Period: It is possible that Co does not 
correspond to a primitive period, that is, it might happen 
that Co can be written in the form: 



UtUUt = U- 



U* 



(6) 



Cn = D' 



(5) 



for some k > 1. In this case, it is not clear to us how 
to define a structural invariance group that leaves this 
property invariant, but the remedy is clear enough: It 
is simpler then to quantize D in the first place. The 
spectral properties of U^ can then be related to those of 

Uco- 

We can now go back to the first two properties, and 
show that a TRI system will have a COE, whereas a 
system with a discrete symmetry will have eigenphases 
thatare divided into statistically independent blocks ac- 
cording to the eigenspaces of Uq. To show the first, we 
point out that the set Q consists of all unitary matrices 
U such that 



Here Ut denotes the anti-unitary map which represents 
time reversal in the quantization under consideration. As 
follows from well-known theorems, it is always possible 
to choose a basis such that Ut corresponds to complex 
conjugation. From (|g) follows that U is a symmetric 
unitary matrix. The action of the quantum equivalent of 
H on il is the following 



Pu:V — y UVU\ 



(7) 



where /7* denotes the transpose of U without complex 
conjugation. Clearly, this action leaves the symmetry 
property of V invariant. Equally, if we act on Uco by 
means of pjj for all U, we obtain as a set, the set of all 
symmetric unitary matrices. We are therefore led to ask 
what measure exists on the set of all symmetric unitary 
matrices such that it is invariant under the action of pu 
given by (R). A standard theorem |13] states that the 
only such measure is the one known as the Circular Or- 
thogonal Ensemble (COE). 

Let us now consider the case of discrete symmetry in 
some detail, following closely the reasoning presented in 
^^- In the case of discrete symmetry, the structural 
invariance group Ti is given by pairs of unitary matri- 
ces acting to the right and to the left of another matrix, 
which is also block diagonal in the same basis. The Haar 
measure of Ti, is clearly the product measure of the Haar 
measures of the various unitary groups restriced to each 
block. From this follows that the ensemble is the prod- 
uct of the various CUE's involved, which is a well-known 
empirical fact, as mentioned in the Introduction. 

Similarly, we had pointed out a difficulty in the case 
where the existence of a discrete symmetry was combined 
with TRI. This can now be treated in the above manner: 
We must determine the structural invariance group Ti 
that respects simultaneously TRI and the symmetry Q. 
In this case, if we describe everything in a basis adapted 
to C, it is not true any more that we can describe T 
through simple complex conjugation. Define Ut to be 
the anti-unitary representation of T in quantum mechan- 
ics. A moment's thought shows that two cases are pos- 
sible: First, Ut leaves invariant the eigenspaces belong- 
ing to the irreducible representations of C . Then, in 
this eigenspace, Ut can be represented as complex con- 
jugation and the ensemble restricted to this eigenspace 
is indeed the COE, using the same kind of arguments 
that lead to the COE when no symmetry is present. 
On the other hand, Ut can act as the combination of 
complex conjugation and the operator interchanging the 
given eigenspace with its time-reversed counterpart. This 
is clearly only possible if some Q in G is not itself TRI, 
that is, if 



TQTQ ^ 1. 



(8) 



The group Ti is then given by the pairs {U, UtU ^Ut)- 
In the block diagonal form, the matrix U has entries 



C/i 

t/2 



(9) 



whereas the matrix UtU ^U' 



T gives 







ul 



(10) 



From this follows readily that the eigenvalues in both 
eigenspaces are degenerate, but that the ensemble for 
each of them is actually the CUE. If we go through the 
above considerations from the point of view of group the- 
ory, we see that this phenomenon can occur whenever 
some irreducible representation of the symmetry group 
is not self-adjoint, that is ,when no basis can be found 
in which all the matrices of the representation are sym- 
metric. The conjugate pair of representations and their 
associated invariant subspaces will then be interchanged 
by time-reversal. 

This was verified numerically in . The specific exam- 
ple considered was a fully chaotic billiard with threefold 
symmetry but without any symmetry axis. The rotation 
by 27r/3 is clearly not TRI and it divides the Hilbert space 
in three invariant subspaces, generated by e*""^ for ra 
equal to one, zero or minus one modulo three respectively. 
Whereas the second is clearly invariant under TRI, the 
other two are interchanged among each other. For these, 
indeed, GUE statistics was clearly observed, whereas in 
the TRI subspace the usual GOE was observed. These 
observations were later confirmed and explained on the 
basis of periodic orbit theory by Keating and Robbins 



IV. FROM EIGENPHASES TO EIGENVALUES 

So far we have only treated the case of a bijective 
canonical map, which can e.g. describe the time evolu- 
tion of a periodically driven system or a Jung scattering 
map [^. An important question, however, is to transfer 
this analysis to flows generated by a time-independent 
Hamiltonian. The immediate problem is that we must 
find a way to account for the fact that the overall density 
of states is not given by that of the corresponding matrix 
ensembles and can in fact be rather arbitrary, whereas 
the fluctuation properties are indeed given by the RMT 
predictions. This difficulty was absent in the earlier cases 
since indeed the density of states is correctly predicted 
to be uniform. 

To do this, we consider the energy-dependent Poincare 
map Ce(p,<i), where the variables refer to phase space 
variables in the Poincare surface of section. The Poincare 
map satisfies the conditions described above: It is bijec- 
tive (up to a set of points that start from the Poincare 
surface but never return to it. By the Poincare recur- 
rence theorem, however, these are of measure zero and 
can therefore be disregarded). Further, the part of phase 
space on which it is defined is a subset of the energy sur- 
face, and as such is automatically compact for a bounded 
system. We therefore deduce that the eigenphases of the 
Poincare map are distributed in the way corresponding 
to the symmetries of the map, which themselves corre- 
spond to the symmetries of the system. It now remains 



to show in which way this distribution of the eigenphases 
carries over to that of the eigenvalues. 

To this end one proceeds as follows: Bogomolny pc| ] 
has shown that a semiclassical quantization condition is 
the following: if £"„ is an eigenvalue, then 



det(l-T(^„))=0. 



(11) 



Thus, if the eigenphases oiT{E) are denoted by exp(i0j {E)), 
then every time a given (j)j{E) goes through zero, E is an 
eigenvalue. It turns out that the whole procedure is only 
semiclassical, as the map T{E) is only unitary in the 
semiclassical limit, but this is not a problem, since this 
limit is in any case the only one we are able to handle. 
Also, outside of the true semiclassical limit, the relation 
between chaos and RMT is much more subtle: in par- 
ticular, one has problems such as (transient) Anderson 
localization, in which chaotic behaviour and randomly 
distributed eigenvalues coincide. One also then has to 
deal with finite tunneling probabilities and other phe- 
nomena associated with the structure of our canonical 
map in the complex plane, which we have left out of con- 
sideration entirely, as our understanding of these is still 
very incomplete. 

It therefore appears that we have reduced the problem 
of determining the spectrum of a Hamiltonian H(j), q) to 
the study of the energy dependence of the eigenphases of 
the quantized version of its Poincare map. Here the tilde 
indicate canonical variables on the complete phase space, 
while canonical variables without tilde are defined on the 
surface of section. To handle the problem, we must first 
know how the Poincare map changes under infinitesimal 
changes of E. To this end, let us consider two nearby 
energies E and E + AE. The map Ce then maps {p, q) 
onto {p', q') and C^+^£; maps (p', q') onto a point (P, Q) 
near to the initial condition. Thus for almost all (p, q) 



P = p + Sp 
Q = q + Sq, 



(12) 



where 6p and Sq are small. The exceptions occur if an 
orbit bifurcates between the two energies E and E + AE. 
Since this transformation is canonical, there must exist 
a "Hamiltonian" Tg (p, q) , such that 



6p = -AE 



dTEip,q) 



Sq = AE 



dq 
drE{p,q) 
dp 



(13) 



that is 



P = p- AE 

oq 

op 



(14) 



Thus the coordinates on the surface of section satisfy 
the Hamilton equations for the canonical map C]^,^^Ce 
where Te{p, q) plays the role of the Hamiltonian and AE 
plays the role of the time. To show that Te{p, q) is the 



time necessary to return to the Poincare surface if one 
starts from (p, q) , we need to use reduced action Se — 
J{pdq — Hdt}. By standard arguments |g^, is easy to 
prove that —Se is the generating function of Ce- Using 
the same arguments, Se+ae is the generating function 
of 



C 



E+AE- 



(15) 



Finally, the differential of the action of the trajectory 
that start on q and ends on Q is given by 



dS = d{SE+AE — Se) 



ah, 



(16) 



(17) 



Finally the time to return to the Poincare surface of sec- 
tion is given by 



TE{ps.qs)^%. 



(18) 



Summarising the results for the Poincare map, the in- 
finitesimal canonical transformation close to the identity 



C 



^eXaeCe = 1 - ^EC-\E + AE))^CiE), (19) 

means the following: consider the "Hamiltonian" TE{p,q)^ 
which is defined as the time necessary to return to the 
surface of section if one starts from (p, q) . This "Hamil- 
tonian" generates a flow on the Poincare surface and 
^E+AE^E is the infinitesimal canonical transformation 
corresponding to following this flow for a "time" /^E. 

If we now follow this through the quantization pro- 
cedure, we obtain the following: Denote the eigenphases 
by <j)j{E) and ^j{E) be the corresponding eigenfunctions. 
Thus 






-.{^^{E)\~^T{E)^\UE)) (20) 



If we now denote by TLe the self-adjoint operator cor- 
responding to 7e, we finally obtain 



i| = w^. 



(22) 



where the overline denotes average over phase space. The 
crucial points to note is that the right-hand side depends 
on a classical energy-dependent quantity. However, it is 
quite independent of the specific wave function and hence 
of J. Since there are N eigenphases on the unit circle and 
they move with a velocity on the order of one, they will 
cross zero at energies which differ by an order of 1/N . 
This is natural, since we have chosen our scale of energies 
to be the classical one. Therefore, from one eigenvalue 
to the next, the velocity at which the eigenphases moves 
hardly changes. Thus, locally at least, the motion of 
the eigenphases is quite rigid. This means that the RMT 
properties of the eigenphases translate directly into corre- 
sponding properties of the eigenvalues of the system. On 
the other hand, it is important to realize that this con- 
stancy does not hold forever. Two effects will eventually 
alter change the velocity: first, the average on the l.h.s. 
of Eq. ( p2| ) will experience a secular change as E changes. 
This corresponds to the secular change in the density of 
states which is usually eliminated by unfolding the spec- 
trum. On the other hand, another effect may well come 
into play even before this secular change becomes notice- 
able. The point is that Eq. ( |2^ ) is only true as a statisti- 
cal statement and there are fluctuations around the mean 
velocity. The most obvious cause for such fluctuations are 
departures of the Husimi distribution from equidistribu- 
tion. Such deviations are well-known to exist, namely 
the so-called "scars" near short periodic orbits. It could 
therefore well be that these accumulated fluctuations ac- 
count for some of the effects due to short periodic orbits. 
To explore this possibility, however, we would presum- 
ably require an understanding of scars which we do not 
have at present. Note also that the long-range stiffness 
found by Berry M is a phenomenon that lies beyond the 
range of validity of the above remarks. Indeed, its on- 
set is at an infinite distance in terms of the mean level 
spacing, as we go to the semiclassical limit. 



V. THE QUANTUM POINCARE MAP ON THE 
BILLIARDS 



dE 



= {^,{E)\nE\^AE))- 



(21) 



Now we must make some key approximations: First, we 
remember that we are in the semiclassical limit, that 
is, that the classical function Te{ps,(Is) is smooth com- 
pared with the ipj{E). This implies that instead of 
using, say, the Wigner distribution in computing the 
l.h.s. of Eq. (8), we can use a smoothed version such 
as the Husimi distribution without great error. Since the 
"tpj (E) are eigenfunctions of a matrix representing a to- 
tally structureless map, their Husimi distributions will be 
spread uniformly all over phase space. This would follow 
from our considerations on structural invariance, but is 
equally confirmed by a rigorous theorem of Shnirelman's 
concerning eigenfunctions of chaotic systems. From this 
one finally gets 



Here we will give a numerical example of the principal 
results of the last section by studying the properties of 
the eigenphases for the rectangle and Bunimovich stta- 
dium with Neumann boundary conditions. The latter is 
defined as the (convex) region enclosed by two semicir- 
cles of radius R connected by two parallel segments of 
length L. This system has been shownJM to be com- 
pletely chaotic. Following Bogomolny pO| and Boasman 
p2[ , we will use the boundary integral method for the 
Helmholtz equation. For a billiard defined by a region A 
with boundary B it is 



-— V2*„(a;)=£;„^„(x) 
Zm 

dn 



X inside A 
X on B, 



(23) 



where ^ is the normal derivative to the boundary 



at point X. Introducing the usual Green function 
Gq{x' ,x;E) for the free Laplacian and using the Green 
identity, on has after integration over the complete re- 
gion, 



1 



^niq') = 



2m 



dq 



dGoiq\q;E) 
dn 



*n(g) - (24) 



d^njq) 
dn 



Go{q\q;E) 



+ {E-En) Goiq\x-E)^n{x)dx, 

J A 

where q and q' denote coordinates on the boundary B. 

Setting E = En and imposing the Neumann boundary 
conditions (other kind of conditions are straightforward) 
the last equation takes the form: 



h" J dGo{q',q;E) 
dn 



^n{q)dq = *„(g')- 



(25) 



Discretising the boundary by N equally spaced points 
{q}i,j=i,...,N and approximating the integral by the trapeze 
rule, the last equation takes the form 

^E^^%^^*"fe)-*"('^0 (26) 



j=i 



dn 



and the eigenvalues of the billiard can be obtained by 
looking for the roots 

det(l-dT,j(£;)) =0 i,j = l,...,N (27) 

when E is varied as parameter. Here d is the distance 
introduced by the discretisation of the integral and 

T.^jiE) = A {n{q) . (q^ _ g,)}i7i(fc|g, - qt). (28) 

The Hankel function is used because in two dimensions 



Go{q.,q,;E) = -^4'^ UJ^lq^^q^lj , (29) 

and n(qi) is the outgoing normal to the boundary at point 
qi. To look for the eigenphases we now solve the eigen- 
value equation for the matrix Tij{E). 

We begin by studying the properties of the eigenphases 
of the quantum Poincare map in a rectangle and in a 
quarter of stadium. The quantum surface of section that 
we used is the boundary of the billiard. The eigenphases 
of Tij{E) for a given value of the energy arc shown in 
Fig. 1. The unitarity of the QPM is better in the high 
energy regime as expected. Also the norms are closer 
to one when the respective phase is zero or tt. This be- 
haviour is shown in Fig. 2 where we plotted the norms 
as function of the phase. 

To study in a detailed way the behaviour of the eigen- 
phases of Fig. 1 as functions of the energy we plot in 
Fig. 3 their phases and norms as function of the energy. 



This figure suggest that the phases for the stadium move 
in a "rigid way" when the energy is varied whereas for 
the rectangle the phases move at least with two velocities. 
To quantify the velocities we plotted in Fig. 4a the dis- 
tribution of velocities for the stadium. From this figure is 
easy to conclude that all eigenphases move basically with 
the same velocity for the stadium. These results shown 
the validity of Eq. ^ (this conjecture was also made by 
Prange [Q) that the eigenvalues in the chaotic case move 
in a rigid way independently of the number of eigenvalue. 
This is not the case for the rectangle (integrable system) 
as we can see in Fig. 4b. 

To study the secular variation of eigenphases we plot in 
Fig. 5 the distributions of the velocities for the stadium 
for 2 different ranges of k. Also, in table I we tabulate 
the mean values as well as the widths of these and other 
velocity distributions. From this it is possible to see that 
these widths are basically constant when k is varied in 
the chaotic case. This is accord with the theoretical pre- 
diction of Doron and Smilansky for a k~^ behaviour of 
the width as function of E. We also show the widths for 
integrable systems in table II. It is seen that, for com- 
parable values of the mean, the widths are significantly 
larger than in the chaotic case. Again, the widths are 
found to be independent of k. The fact that the widts 
are larger in the integrable case is therefore quite robust 
under variations of the energy. 

The numerical results of this section thus confirm the 
theoretical predictions: for the completely chaotic case 
the eigenphases move almost with the same velocity and 
thus its statistical properties translate to those of the 
energies. Nevertheless this conclusion does not hold for- 
ever. When the eigenvalues have given a complete turn 
on the unit circle, large correlations must occur. These 
long range correlations are well known: the saturation 
in the spectral rigidity (equivalently in the Dyson-Mehta 
statistics) related to the shortest periodic orbit [0, and 
the periodicity in the Fourier transform of the two levels 
cluster function, b2{t) for Rydberg molecules pi]. Also 
eigenphases with small velocities have been found in the 
stadium. These correspond to the bouncing ball states. 
Note finally that the techniques proposed by Lombardi 
and coworkers [p6| for Rydberg molecules in the frame- 
work of Multichannel Quantum Defect Theory can be 
interpreted from the same point of view and provide an 
exactly unitary Quantum Poincare map. These maps dis- 
paly similar properties as those described above for the 
Poincare map of the stadium. In particular, the stiffness 
in the motion of the eigenphases for the case of chaotic 
motion is also observed in this system [pSi . 



VI. CONCLUSIONS 

We have presented a systematic approach to the con- 
nection between RMT and individual dynamical systems. 
This connection is in a sense of a probabilistic type: it 
rests basically on conclusions of the form: let a given sys- 
tem be a typical representative of a certain ensemble; If 
the elements of this ensemble have a given property with 



probability one, the original system presumably also has 
the stated property. While this line of argument is wide 
open to criticisms from the mathematical side, it is unde- 
niably useful at the heuristic level. Further, similar lines 
of reasoning are frequently used in physics — for example 
in statistical mechanics — when trying to apply an ensem- 
ble description to an individual system. A more genuine 
concern concerns the construction of the ensembles: as 
we are not able to construct ensembles on the set of all 
canonical transformations, we must first construct a set 
of canonical transformations at the purely classical level 
and then translate this into quantum mechanics in order 
to obtain a reasonable candidate for a measure. On the 
other hand, it may be possible — and we would like to sug- 
gest that this would be important — to find an invariant 
measure on the set of all bijective canonical transforma- 
tions on a compact phase space. If this were possible, we 
could define natural classical equivalents to the CUE and 
the COE. One could then argue in an entirely classical 
way that a given map with chaotic dynamics is a typical 
element of such a classical ensemble of canonical maps. 
Its quantization would therefore belong to the associated 
quantum ensemble, whence the required spectral proper- 
ties follow. 

Another limitation is the purely semiclassical nature 
of the analysis. It does not, for example, apply to quan- 
tum localization: as long as the phase space T is com- 
pact, localization is only transient, and therefore outside 
the immediate range of application of semiclassics. Non- 
compact phase spaces, on the other hand, also present 
problems relating to the existence of an invariant mea- 
sure, not only in the classical, but also in the quantum 
case. 

On the other hand the translation of the fluctuation 
properties from the eigenphases to the energy spectrum is 
particularly transparent. Thus the fluctuation properties 
of the energies are the same as those of the Gaussian en- 
sembles with probability one. This approach shows that 
not only do the RMT predictions hold for the two-point 
function at intermediate energy distances, but that they 
should hold at all energy scales and for all correlation 
functions. Also the saturation of the long rage stiffness 
for the energy spectrum is understood in terms of com- 
plete turns of the eigenphases on the unit circle. Finally 
we want to mention that the properties of the velocities of 
the eigenphases promise to be a new quantum signature 
of chaos. 
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FIG. 1. Eigenvalues of the operator associated to the quan- 
tum Poincare map for a quarter of stadium of size L = 20 y 
R = 40. We have taken iV = 58 and k = 0.49. Note that the 
unitarity is better for phases close to and n. 
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FIG. 2. Norms of the eigenvalues of Fig. 1 as a function of 
their phases. We plotted only those with a norm larger than 
0.8. 
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FIG. 3. a) Norms and b) phases of the eigenphases of the 
quantum Poincare map for the same stadium for A*' = 58 as 
a function of k. 
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FIG. 4. Velocity distribution for a) the same quarter of 
stadium for k ~ 0.45, A'' — 290 and b) the rectangle with an 
aspect ratio given by the golden mean. We use A*' — 300. We 
only used the eigenphases whose norm is greater than 0.6. 
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Center 


Width 


0.02 


515.438 


77.142 


0.03 


516.801 


85.977 


0.04 


519.923 


82.403 


0.05 


519.559 


78.797 


0.06 
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73.579 


0.07 


521.260 


72.360 


0.08 


522.025 


72.095 


0.09 


522.057 


71.590 



TABLE II. Centers and widths of the velocity distributions 
of the same rectangle for various values of k and N = 300. 
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FIG. 5. Velocity distributions for the same quarter of sta- 
dium with iV = 290 and a) fc ~ 0.02 and b) fc ~ 0.09. 
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42.575 
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33.713 



TABLE I. Centers and widths of the velocity distributions 
of the same quarter of stadium for various values of k and 
N = 290. 
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